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Abstract 

We consider the first order periodic systems perturbed by a 2N x 2N matrix-valued 
periodic potential on the real line. The spectrum of this operator is absolutely con- 
tinuous and consists of intervals separated by gaps. We define the Lyapunov function, 
which is analytic on an associated N-sheeted Riemann surface. On each sheet the Lya- 
punov function has the standard properties of the Lyapunov function for the scalar 
case. The Lyapunov function has branch points, which we call resonances. We prove 
the existence of real or complex resonances. We determine the asymptotics of the pe- 
riodic, anti-periodic spectrum and of the resonances at high energy (in terms of the 
Fourier coefficients of the potential). We show that there exist two types of gaps: i) 
stable gaps, i.e., the endpoints are periodic and anti-periodic eigenvalues, ii) unstable 
(resonance) gaps, i.e., the endpoints are resonances (real branch points). Moreover, we 
determine various new trace formulae for potentials and the Lyapunov exponent. 

Keywords: periodic systems, spectrum, high energy asymptotics. 



1 Introduction and main results 

Consider the self-adjoint operator K, acting on the space L 2 (R) 2N and given by 

Ky = -u iv > + v tv , j^C" _y , 

here and below we use the notation (') = d/dt and Ijv is the identity N x N matrix; v is the 
complex 1-periodic N x N matrix and V = V t belongs to the Hilbert space Jf? given by 

^ = { V = (v* o) ' v = vT= to*W}&=i' 1 e R / Z > ll^ll 2 =f Q Tl V ? dt < °°}- 

*Institut fur Mathematik, Humboldt Universitat zu Berlin, Rudower Chaussee 25, 12489, Berlin, Ger- 
many, e-mail: evgeny@math.hu-bcrlin.de 



1 



Note that Rev and Imti are self- adjoint. Without loss of generality we assume 

y=[ V t 2 dt = %®%, % = diag{v u ...,v N }, < v x < u 2 ^ ... ^ u N , (1.1) 
Jo 

here v\, .., u N are the eigenvalues of f, see Sect. 2 for the proof. 

It is well known (see [DS] p. 1486-1494, [Ge]) that the spectrum o(K) of K, is absolutely 
continuous and consists of non-degenerated intervals a n , n G Z. These intervals are separated 
by the gaps <7 n = (z~,z£) with the length |g n | > 0, < n < iV+, where — oo ^ iV~ < 
iV+ ^ oo and N g = — N~ — 1 is a total number of the gaps. 

Introduce the fundamental 2N x 2iV-matrix solutions ip(t,z) of the equation 

- U-rf + V t i> = zi>, z G C, ^(0, *) = I 2 n, 

and the monodromy 2N x 2iV-matrix ^(l,z). The matrix valued function ^(1, •) is entire. 
An eigenvalue of ip(l,z) is called a multiplier of /C: to each of them corresponds a solution 
/ of —%J\f + 14/ = zf with /(£ + 1) = T(z)f(t),t G [0, 1). They are roots of the algebraic 
equation D{r,z) = 0,r, z G C, where D(t,z) = det(ip(l,z) — t/jv), r, 2; G C. Zeros of the 
function det(.M(2) — I N ) (or det(.M(.z) — ijv)) are periodic (or anti-periodic ) eigenvalues. 

There exist many papers about the first order periodic systems iV > 2: Gel'fand and 
Lidskii [GL],Gohbert and Krein [GK], Krein [Kr], Potapov [Po], [YS]) and we mention new 
papers of Gesztesy and coauthors [CG],[CHGL], [GKM]. The basic results for direct spectral 
theory for the matrix case were obtained by Lyapunov and Poincare (see [GL],[Kr],[YS]). 
Theorem (Lyapunov, Poincare). For each (V, z) G Jtf? x C the matrix-function ^>(l,z) 
satisfies: 

r\i,-) = -J4> T (i,-)J, (1.3) 

D(t,-) =t 2N D(t-\-), any r ^ 0, where D(t, •) = det(^(l, •) - rI 2N )- (1.4) 

a(K) = {z G C : \t(z)\ = 1 for some multiplier r(z) of ^(M)}- (1-5) 

If for some z <E C (orzEM) t(z) is a multiplier of multiplicity d ^ 1 , thenr~ x {z) (orr(z)) 
is a multiplier of multiplicity d. Moreover, each ip(l,z),z G C, has exactly 2N multipliers 
T^ zl (z),j = 1, ..,N. If t(z) is a simple multiplier and \r(z)\ = 1, then r'(z) 7^ 0. 

The eigenvalues of ^(1, z) are the zeros of the equation D(t, z) = 0. This is an algebraic 
equation in r of degree 2N with coefficients, which are entire in z G C. It is well known (see 
e.g. [Fo]) that the zeros Tj(z),j = 1,..,2N of D(t,z) = are (some branches of) analytic 
functions of z with only algebraic singularities: the zeros Tj(z),j = 1, ..,27V constitute one 
or several branches of one or several analytic functions that have only algebraic singularities 
in C. Thus the number of eigenvalues of z) is a constant N e with the exception of some 
special values of z (see below the definition of a resonance). In general, there is a infinite 
number of such points on the plane. If the functions Tj(z),j = 1, .., 2iV are all distinct, then 
N e = 2N. If some of them are identical, then we get N e < 2N and z) is permanently 
degenerate. 

Introduce the matrix-valued function Sf(z) = j(ip(l,z) + ^ -1 (1, z)), z G C and the 
function $(z,v) = det (Jf(z) — uI 2 n),z,u G C. Each zero of $(v,z) has multiplicity ^ 2 
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and define the Lyapunov function by Aj(z) = \{rj{z) + r^ 1 (z)), j = 1, ..,N. The Riemann 
surface for the multipliers Tj(z),j G = {1, .., N} has 2N sheets, see ( 11. 4ft . If N — 1, then 
it has 2 sheets and the Lyapunov function is entire. Similarly, in the case N ^ 2 it is more 
convenient for us to construct the Riemann surface for the Lyapunov function, which has N 
sheets. We need the following results from [K4]. 

Theorem 1.1. Let V G M' . Then there exists an analytic function A s , s — 1, .., Nq ^ N on 

the N s -sheeted Riemann surface £% S) N s ^ 1 having the following properties: 

i) There exist disjoint subsets u s ,s = 1, .., N Q ,[ju s = Nn such that all branches of A s ,s = 
1, 2, .., No has the form Aj(z) = \{jj(z) + r~ (z)), j G u s . Moreover, for any z, r G C the 
following relations hold true: 

N 

det(Sf(z)-vI 2N ) = ]l& l 8 (v,z), $ a (v,z)=Y[(u-A j (z)), z,ueC, (1.6) 

Aj(z) = cosz + o(e |Im " 1 ) as \z\ -»• oo, (1.7) 

where the functions $ s (i/,z) are entire with respect to u, z G C. Moreover, if A, = Aj for 
some i G u>k,j G oj s , then = $ s and Ak — A s . 

ii) (The monotonicity property). Let some Aj,j = 1, ..,N, be real analytic on some interval 
Y = (at, (3) C R and -1 < Aj(z) < 1 /or any z G Y\ T/ien A^(^) 7^ /or eac/i z G F. 

The functions p,p s given by (11.81) are entire, 

No 

p=i[ps, p.(o= n (^(o - A i(-)) 2 - 

iv) The following identity holds true 

a(K)=R\uf =1 g n , g n = (z~ , z+), N g < n < N+ (1.9) 

where each gap g n = (z~, z£) is a bounded interval and z^ are either periodic (anti-periodic) 
eigenvalues or real branch points of Aj (for some j = 1, .., N) which are zero of p (below we 
call such point a resonance). 

Remark. 1) In the case of 2 x 2 system the monodromy matrix has exactly 2 eigenvalues 
r, r _1 . The Lyapunov function |(r + r _1 ) is an entire function of the spectral parameter. 
It defines the band-gap structure of the spectrum. By Theorem 1.1, the Lyapunov function 
for 2N x 2iV-matrix operator /C also defines the band-gap structure of the spectrum, but it 
is the N-sheeted analytic function. 

2) We have the following asymptotics (see Sect. 3) 

si n z / 'pi ^ m z I \ 
Aj{z) = cos z + -^Z Vj + 0(—^- ), ifV'e^ (1.10) 

as \z\ — > 00, j G Nat. Then firstly, p is not a polynomial since p is bounded on R. Secondly, 
if Vji 7^ Uj,j' 7^ j, then (ll.lOp implies Aj/ 7^ Aj. 
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3) In the case N = 2we determine Aj, A 2 , p in terms of the traces of the monodromy matrix. 
Using (TOD, we have D(r, •) = r 4 - T x t 3 + §(T 2 - T 2 )r 2 - 2Vr + 1, which yields 

D(r,.)=(r 2 -2A 1 r+l)(r 2 -2A 2 r + l), A, = | + ^, A 2 = | - ^, (1.11) 

see [BBK], where p = — ^- and T m = Trip(m, z), m — 1, 2. 

Definition. v4 zero Zo £ C o/ p gwen fry (11.81) zs a resonance o//C. 

The main goal of this paper is to describe the spectrum of /C and to determine the 
asymptotics of gaps and resonances, periodic and anti-periodic eigenvalues at high energy. 
We show that all resonances are real at high energy. Moreover, we prove the existence of 
complex resonances for some specific periodic potential. We have to underline that in the 
case of large N the resonances create the gaps in the spectrum of periodic operators, see 
Theorem 1 1 . 2 1 and remark after Theorem II .31 If N = 1, then 2-periodic eigenvalues create the 
gaps in the spectrum. In the present paper we use some techniques from [K4] and [BBK], 
[BK], [CK]. 

The periodic eigenvalues (n is even) satisfy 

.. < z N 2 ' + sc z G { + <: z Q 2 + <: ... ^ 4 + < Z Y < z l' + < ••• < z n~ < z n + < z t'~ < > ( L12 ) 
* , ' * v ' 

n=0 n=2 

the anti-periodic eigenvalues (n is odd) satisfy 

.. sC z N h+ ^ z\~ sC z{' + s= ... ^ z]f ^ z^ + ^ zf~ ^ z 3 { + ... ^ z% + zl~ ^ ... (1.13) 

V v ' V v ' 

n=l n=3 

and they have asymptotics 

z"* = im + o(l) as n^±oo, j G N N = {1, 2, .., N}. (1.14) 

If V = 0, then these eigenvalues have the form z n ^ = irn, (n,j) 6Zx Nat. 
Let = J Q Vle i27rntJl dtj We formulate our first main result 

Theorem 1.2. Let V, V G M 1 and let CT^ij G be eigenvalues of the matrix 'f — iJ\V^. 
Then the periodic and anti-periodic eigenvalues have the following asymptotics: 

z»' ± = 7m +k_ + o(n- 2 ), j en N as n^±oo. (1.15) 

Assume that i/j ^ v.y for all j ^ j' G uo s for some s = 1, .., Nq. Then the function p s has the 
zeros z^jCn = < j',j,f G w s , n G Z ; which are real at large \n\ and satisfy 

^± = vm + ^±^ + o(^ + ^-Y a = {j,j') as n^±oo. (1.16) 
47m V n n z / 
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Let in addition U\ < .. < u^. Then for each s — 1, .., -/V and /or /arge n — > ±00 iaere 
esiste a system of real intervals (gaps) g™ = (z"~,z™ + ) such that 

= « = & A ^" e < < < < - < ^ < 

(-l) n Aj(z) > l,z G and A» = A f (z),z G if j" (1-17) 

£aca branch Aj is real and is analytic on the set (7m — |,7m + |) \ Up^yg^j and is not 
real on U p ^-g™j. 

iij If z™~ 7^ z™ + /or some a = 7^ j', i/ien zs i/ie simple branch point (resonance) 

for the functions Aj, A^. If z™~ = z™ + , then Aj, Aj/ are analytic at z^ . 
Hi) The following asymptotics hold true: 

zt = vrn + ^ + ^ Kal + O (fc| + i) , v^ a =j\' a {t)e~^ nt dt, a = (j,f). (1.18) 

Remark. 1) iV 9 = iV+ — iV~ — 1 is the total number of gaps in the spectrum of /C. 
2) If //!<..< i/jv; then there exist infinite number of resonances , a = 7^ f which 

form the gaps in the spectrum of /C, see below Theorem 11.31 Roughly speaking, resonances 
"form" the gaps, the number of periodic and anti-periodic eigenvalues is less than the number 
of resonances. Thus there exists big difference between N — 1 and large N. In the first case 
the endpoints of the gaps are 2-periodic eigenvalues. In the second case, roughly speaking, 
the endpoints of the gaps are resonances. 

In the second main result we describe finite band potentials. 

Theorem 1.3. Let V, V G and let u x < ... < u N . 

(i) If the identity u% + u^ = u 2 + u^-i = ■•• = v n + ^1 is not fulfilled, then \Ng\ < 00. 
(%%) If v\ + z//v = ... = vn + vi holds true and there exists a sequence of indicies rik — > ±00 
such that \v' na \ 2 + (rifc | 1 = o(|t>^J) as k —> ±00, for each a = (j, N + 1 — j),j G Nn, then 
N± = ±00. 

Remark. 1) Consider v = diag{t>i, V2, ..,vn}, i.e., the case when "variables are separated". 
The transformation %:y= ..,y 2 N) T -> (2/1, 3/iV+i> 2/2, 2/JV+2, 2/27v) T gives 

%K%* = ©f JC jf Ki = -ihj t + Vj, h = Q ^ ) V 3 = Q . 

The operator JCj for the case N = 1 is well studied [YS], [Kl-3], [Mi]. We have Uj = 
Jo \ v j(t)\ 2 dt ^ 0. If Uj < Up for some j < p, then the number of gaps is N g < 00. 

2) Note that the condition |i4a| 2 ~H n l -1 = °(l u nal)j a = C7>-W + l ~ j)->3 ^iv asn-> ±00, 
holds true for "generic" potentials V, V G J^. This yields the existence of the real resonance 
gaps (z%~, z™ + ) at high energy. The coefficients v' na , a = (j, iV + 1 — j), j G Nat (the second 
diagonal of the matrix v) "create" the gaps. 

Example of complex resonances. Consider the operator AT^ T = —iJi^ + V u ^ T ,u = 
1, |, |, .., r G R acting in L 2 (IR) 4 , where the real periodic potential V^ jT is given by 

^ = («" "o r ) • = - ( X t ) • s e R+ ^ N - K e C(T) - (L19) 
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pi pm pm 00 1 

/ \b v (t)\dt=\, / b v (t)f(t)dt^ 5 per (t)f(t)dt, ^)=r^-n--), (1.20) 
Jo Jo Jo z 

as z/ — > 0, for any / G C(0, m),m G N. 

If r = 0, then the operator /C^o = — zJi4 + K,o has the constant potential V v $. In 
this case there are no gaps in the spectrum of K. v> o and all resonances are given by r° = 
Tin + ■f^r,n G Z \ {0} with multiplicity 2. We show that there exist the non-degenerated 
resonance gaps for small t, v. In this example some resonances are real and some are complex. 

Proposition 1.4. Let a potential V T}1/ satisfy ( |1.19p . (I1.20p . Then for each large integer 
n ^ 1 + a there exist sufficiently small u, e > such that the following statements hold true: 
i) Each function p(z, V TjV ), r G (—£, e) in the disk z G D no = {|z| < 7m Q + 1} has exactly 4n 
simple zeros r^ u (r),l ±n ^ no, where r^ u (r) is analytic function in the disk r G {|t| < e} 
and r„j,(0) = r° and the following estimates hold 

ri„(r)=r°±r(V^: + o(l)) « r - 0, ^ere | £j < | , (1.21) 

/or some constant R n and if r G (— e, e), |n| > j- ; i/ien (r~ J/ (r), ^^(r)) Ctisa gap. 

Each function D(l, z, V T<U )D{— 1, z, V^„) m £/ie dzs/c IED no nas e2;ac%4n +4 ^eros ^ m (r, z/), 
— no ^ n ^ no, m — 1, 2 and ine asymptotics z^ m (r, v) = z®'^ + o(l) ao/ds as r — > 0. 

Remark. 1) If < a < 27r, then V^) has only real roots r^(r) in each large disk D no 
for sufficiently small t,u. If 1 < j-, then V^j,) has at least two non-real roots rf^(r) 
for small r, z/. 2) We show that operator )C Ttll has new gaps (so-called resonance gaps). The 
endpoints of the resonance gap are the branch points of the Lyapunov function, and, in 
general, they are not the periodic (or anti-periodic) eigenvalues. These endpoints are not 
stable. If they are real (see (ll.2ip ). then we have a gap. If they are complex (0 < n ^ j-), 
then we have not a gap, we have only the branch points of the Lyapunov function in the 
complex plane. 3) We have a similar complicated distribution of other resonances, which are 
poles of S-matrix for scattering for Schrodinger operator with compactly supported potentials 
on the real line see [K5], [Z]. 

We consider the conformal mapping associated with the operator /C. Introduce 
the simple conformal mapping 77 : C \ [—1, 1] — > {( G C : \(\ > 1} by 

r)( z ) = z + Vz 2 - 1, zeC\[-l,l], and n(z) = 2z + o(l) as \z\ -> 00. (1.22) 

Note that n(z) = rj(z), z G C \ [—1, 1] since T](z) > 1 for any z > 1. Due to the properties 
of the Lyapunov functions we have |r/(A s (^))| > 1,£ G = {( G & s : Im£ > 0}. Thus 
we can introduce the quasimomentum kj,j = 1,2,..,N (we fix some branch of arccos and 
Aj(z)) and the function qj by 

kj(z) = arccos Aj (z) = ilogr)(Aj(z)), Qj{z) = lmkj(z) = log \r)(Aj(z))\, (1-23) 

z G Mq = C + \(3 + , j3 + = U/3eB A nc+ \Pi @ + ^°°) where £>a is the set of all branch points of the 
function A. The branch points of k m belong to £>a- Define the averaged quasimomentum 
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k, the density p and the Lyapunov exponent q by 

N 

,\ 



1 N 

k(z)=p(z)+zq(z) = -J2k 3 (z), q(z)=lmk(z), z G . (1.24) 



i 



For the function k(z) = p(z) + iq(z), z G C + we introduce formally integrals 

Qn = - I t n q(t)dt, I S n = - [ t n q(t)dp(t), % = - [[ \k{ n) (z)\ 2 dxd yj (1.25) 

n — 0, 1, 2, here and below fc(o) = k — z, km = zkm, z = x + iy G C. Let C us denote the class 
of all real upper semi-continuous functions h : R — ► R. With any h G C us we associate the 
"upper" domain K(h) = {k = p + iq E C : q > h(p),p G R}. We formulate our last result. 



iC 



Theorem 1.5. i) Let V G Jtf. Then the averaged quasimomentum k = jt J^i % is analyt 
in C + and k : C + — >■ k(C+) = K(/t) zs a conformal mapping for some h G C us . Furthermore, 
for some branches kj,j = 1.., iV the following asymptotics, identities and estimates hold true: 

k(z) — z = — — _ as ^1 — > oo j if y > r\x\, for any r > 0, (1-26) 

Qo = / D + ^ = ^T, (1-27) 

g| CT(JV) = 0, < g 2 U (1) u 9 < 2Q , wfcere ct(at) = {^ G R : Ai(z), .., Ajv(z) G [—1, 1]}, 
<T(i) = {z G R : A m (z) G [-1,1], A p (z) £ [-1,1] some m,p = 1,..,N}, 

n 

Let aditionly V G Jf 7 . JTien i/ie following asymptotics, identities hold true: 

U ( \ Qo Ql Q2 + 0(l) . . . , 

k(z) — z = as \z\ — > oo, if y > r\x\, for any r > (J, (1.29) 

z z 2 z 6 

Qi = wf Q ~ iJ ^ v ^ = I ? +I *-^ L = m I! ( K/2 + v ') dL (L30) 

A priori estimates for various parameters of the Dirac operator (the norm of a periodic 
potential, effective masses, gap lengths, height of slits, action variables for NLS and so on) 
were obtained in [KK1-2], [Kl-3], [Mi] only for the case N — 1. In order to get the required 
estimates the authors of [KK 1-2], [Kl-3], [Mi],... used the global quasi-momentum as the 
conformal mapping, which was introduced into the spectral theory of the Hill operator by 
Marchenko-Ostrovski [MO] . 

The mapping k : C + — > K(/i) is illustrated in Figured] for the case N = 2. The integral 
Iq ^ is the area between the boundary of K(/i) and the real line. In Figure 1 the upper 
picture is a domain K(/i) and the points A = k(A),B = k(B),... . The spectral interval 
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(A, B) (with multiplicity 2) of the z-domain is mapped on the curve (A, B) of the k-domain, 
the interval (a gap) (B, C) of the z-domain is mapped on the vertical slits, which lies on the 
line Refc = 0. The spectral interval (C, D) (with multiplicity 2) of the z-domain is mapped 
on the curve (C, D) of the k-domain. The spectral interval (D, E) (with multiplicity 4) of 
the z-domain is mapped on the interval (D, C) of the k-domain. The case of the interval 
(E, J) is similar. The resonace gap (K, L) of the z-domain is mapped on the vertical slits, 
which lies on the line Re A; = 2n. In fact we have the graph of the function h(p),p G R, 
which coinsides with the boundary of K(/t). 



fe-plane 
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Figure 1: The domain K(/t) = k(C + ) and the graph of the function h 

We describe the plan of our paper. In Sect. 2 we obtain the basic properties of funda- 
mental solution ip(t, z). In Sect. 3 we determine the asymptotics of the fundamental solution 
ij;(t, z) for the case V, V G Jff. In Sect. 4 we determine the asymptotics of the Lyapunov 
function and multipliers at high energy and prove Theorem II. 2[ 11.51 In Sect. 5 we prove 
Proposition 11.41 In Sect. 5 we determine the asymptotics of detJf(z) as Imz oo. 

2 The fundamental solutions 

In this section we study ip. We begin with some notational convention. A vector h = 
{h n }f G has the Euclidean norm \h\ 2 = \ h n \ 2 , while a. N x N matrix A has the 
operator norm given by \A\ = supm =1 \ Ah\. Note that \A\ 2 ^ Tr A* A. 
Recall the identity: if A is iV x iV matrix, then 

det(A-vI N ) = (-l) N J2w N ~^ oo = l, a 1 = -A 1 , 2 = - M ' A ' ^\ (2.1) 

o 

and .., 4>j = -jEiAfli-*) ••, a>N = detA, A m {z) = TiA m . 
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Below we need the identity 

J 1 V = -VJ 1 , e zJl V = Ve~ zJ \ (2.2) 
for any z G C. The solution of the equation ( 11.21) satisfies the integral equation 

if,(t lZ ) =tfo(t,z) -i [ e izJl{t - s) J 1 V s ^{s,z)ds, Mt,z) = e iztJ \ t > 0, z G C. (2.3) 
Jo 

It is clear that Eq. (I2.3P has a solution as a power series in V given by 

4>(t,z) = Y>n(M), ^n(M) = ~i I e izJ ^J 1 V s ^ 1 (s,z)ds, n^l. (2.4) 

Using ([2J]), (J2I3D we have 

if> 1 (t,z) = -i [ e lz - h{t - s) J 1 V s e izs - h ds = -i [ e lzJl{t ~ 2s) J^ds, (2.5) 



ih = -i f e izJ ^ t -^J l V tl Mti,z)dt 1 = [ dh f 1 e lzJ ^~ 2t ^ 2t W tl V t2 dt 2 . (2.6) 
Proceeding by induction, we obtain 

pt ptln-\ 

^ 2n (t, z)= dh... e izMt-2t 1+ 2t 2 ...+2t 2n ) Vti _ _ _ Vt2 jt 2nj (2.7) 



^2n+l(t, Z) = ~l dh... e izJi(t-t 1+ 2t 2 ...-2t 2n+1 )j iVti Vt2n+1 dt 2n+1 . (2.8) 

Jo Jo 
We need the following results from [K4]. 

Lemma 2.1. Let V G Jtf. For each z G C t/jere exists a unique solution ip of Eq. (12.31) given 
by ( 12.41) and series (" 12.41) converge uniformly on bounded subsets of R x C x . For eac/i 
t ^ the function if)(t, z) is entire on C. Moreover, for any n ^ and (t, 2) G [0, 00) x C 
t/te following estimates and asymptotics hold true: 

JImz|t , pt . n 

< H / lV.|ds) , (2.9) 



711 

n— 1 

n! 

D iztJi _ / t| Im^h 



wt,z)-r^,(M)i < (^iyi)" e tiim^/> s |d S) (210) 



V>(t,z) -e i2tJl =o(e t|lm21 ) as |z| -> 00, (2.11) 

uniformly on bounded i 6 R. If the sequence V v — > V weakly in Jf?, as v — > 00, then 
tp(t, z, V u ) — > tp(t, z, V) uniformly on bounded subsets o/f xC. 
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Below we need the simple properties of matrices a,b,c E A = |a = f ® a ^ ^ : 
ai, a 2 is N x N matrix| given by 

abc, J x a, e zJl aeA, all z E C, (2.12) 

ah = (— axh + a 2 b 2 ) + (axb 2 - a 2 &i)Ji, ( 2 -13) 

Tr a = 0, Tr Ja n = 0, n ^ 0. (2.14) 
For any matrices A, B the following identities hold 



Tr AB = Tr BA, Ti A = Ti A*. (2.15) 
Using (12.6p - fj2.8p we define the function 

Tj )X (z, V) = Tr f V S1 [ 1 V S2 e izJl ^- 2si+2s2) ds, (2.16) 
Jo Jo 

T j , n (z,V) = Tnknti,z) = r ft /•••/ e izJ ^- 2s ^ 2s2 - +2s ^V sl ...V S2n ds, n > 2, (2.17) 

Jo Jo 



where s = (s X) .., s 2 „) G to2 " 



Lemma 2.2. Lei V G M 3 . The functions Tj(-,V),j = 1,2,..,N are entire on C and 
T(z, V) G R /or a// zGl. Moreover, the function T(z, V) satisfies 

T j {z,V) = T j {z,-V) = 2Ncosjz + J2 T j,n(z, v ), ( 2 - 18 ) 

IT,-^,^)! ^ 2iVe j(|Im * l+l|y|l) , (2.19) 

Tj(z,V) = 2Ncosjz + o(e jllmzl ) as \z\ -> oo, (2.20) 

and Tr ^n+i (^j 2) = 0- Series (12.181) converge uniformly on bounded subsets o/Cx J^ 5 . //a 
sequence V v converges weakly to V in J4? as v — > oo ; i/ien Tj(zA^ u ) — > Tj(z,V) uniformly 
on bounded subsets of C. 

Proof. By Lemma I2.2[ series (I2.18P converge uniformly and absolutely on bounded subsets 
of C x TC. Each term in (12.181) is an entire function of z, then Tj is an entire function of z, V. 
Moreover, if the sequence V v converges weakly to V in Ti, as v — > 00, then Tj(z,V u ) — > 
Tj(z, V) uniformly on bounded subsets of C. 

We have 7} = Tri]){j,z) = Tr J2n>o i'nU, z ) an d Trip (j,z) = 2N cos jz for j ^ 1. The 
estimate fl23]) yields (l2T9|) and (Tjmf gives (l2^20|) . 

The first relation in (12.121) yields V tl . . . V t2n+1 G .4 for any t 1; . . . t2 n +i G R. Then relations 
in (I2TT2"]) give Tr e izJl(t - 2tl+2t2 -- 2t2n+l) JxV h .V t2n+1 = 0, for any t, t a , . . . t 2 n+i e R, which 
together with (12.111) implies Tr ■02n+i (^ z) — 0. 

Below we will show that T(z, V) EM. for all z E R (see ( 12~28D ). ■ 
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We will show (II. ip . Consider the self-adjoint operator Jy' + Q(t)y,J 
acting on the Hilbert space (B 2N L 2 (M>), where the real matrix Q is given by 



I N 
-In 



fix tt 2 



^=(^J^ ): ^2 = ^2, ^i = ^l, \\n\\ 2 =j TrVl 2 (t)dt < oo. (2.21) 

Let M(t, z),t eRhe the fundamental solution of the equation JM' + VLM = zM, M(0, z) = 
I 2N - Note that M(t, z) is real for zGl 

Define the unitary matrix — -^(-A + iJ) = % 2 = I 2 n- Using the identities 

J 2 = ( J ^ j , <^J<^ = _iJi, ^ Ji^ = iJ, <^J 2 ^ = ~J 2 , = Ji^i+J2^2, (2.22) 

we deduce that M c = satisfy the equation -U X M' C + Q C M C = zM c , M c (0, z) = I 2N , 

where Q c is given by 



1 r l r l 

2/ 



Thus we obtain 

Q 2 c (t)dt = f{ Q) %, %= I co{t)cu*(t)dt, %= I u*(t)u{t)dt, (2.24) 
o Jo Jo 

and 

% = E%E* = y 2 ^0, % = diagK .., u N }, (2.25) 

for some unitary matrix E and the diagonal matrix Define the unitary matrix S = E@E. 
The function ip(t, z) = <S*M c {t, z)S satisfies —iJ\ip' + Vip = zip, ip(0, z) = I 2 n, where 

<g = E®E ) V = £*n c £= ( JJ V ) , v = E*coE, (2.26) 



i 

2 



y^dt = % © r , (2.27) 

Tr^(t,z) = TrM(t,z), (2.28) 

which gives V G and Tr-^(1, z) G R for all zGi 

It is well known that for real Q we have M(t, z)JM(t, z) T = J (see [GL], [YS]). Then 
M c = WMfy and fl2T22|) give 

-U x = M c (l, z)(-iJi)W& T (M D (l, z) T )W T W, (-iJ^W^ = -J, 

which yields 

M c (l, z)JM c (l, z) T = J. (2.29) 
The similar arguments and ip = <§*M C § imply 

^g*Jg* T iP T g T = g*J, g*Jg* T = J, 
which yields ipJip T = J and (11.31) is proved for V G Jtf. 



11 



3 Estimates of ip for the case V E Jrif 

Sect. 2 does not give the needed estimates of the fundamental solution ip at high energy. In 
order to determine the asymptotics of ip we will do some modification. Define the integral 
operator K and the matrix- valued function a t (z) by 

(Kf)(t) = Je t . s WJ(s)ds, W = -i.hV 2 - V, e t = e t (z) = e UzJ \ a t (z) = I - ^. (3.1) 

where e t = e t (z) = e ltzJl . Introduce 

V~'(z) = [ Vje t2tzJl dt, V*(z) = [ V t 3 e i2tzJl dt, V= I V t 2 dt. (3.2) 
Jo Jo Jo 

Lemma 3.1. For each (z, V) G (C \ {0}) x 7i the solution ip = a~ 1 ^ao, where ^ satisfies 
* = + eKa' 1 ^, 1 & = aipa 1 , a = I - eV , e = — , (3.3) 

* = + Yl *»> = ^(^«~ 1 )>o, (3.4) 

where series (13.41) converge uniformly on bounded subsets ofM. x (C \ {0}) x 7Y. Moreover, 
i/sup tgR |Vt| ^ £aen /or any j — 1, m G N iae following estimates are fulfilled: 

|* B (f,z)|<__ / |W s |d S , (3.5) 

P izJ i / ^ 1 / — - f 1 \\ JImz| 

$!(!,.) = -— ( ! / 1 f + y'( 2 ) + -( l J 1 ^(,) + jf ^y tC ftjj+o(_ ), (3.7) 



* 2 (M) = -j^J dtjy^-tJ^V^-tVle^V^-VlenV^ys + Oi^) (3. 



as \z\ — > oo and where et = e lztJl . 

Proof. Using (12. 21) . H 2 . 3 f) . e = j- and integrating by parts we get 



V>(t, «) - e tztJl = -i j J 1 e izJ ^ t - s) V s i){s ) z)ds = -i J e tzJl ^ 2s) J x \V s e- izsJl ^{s, z))ds 

= ee izMt ~ 2s) (v s e- lzsJl ^{s,z)^ * - e J e LzJ{t ~ 2s) (v s e~ izsJl ^{s, z)^j ds 
= e(v t ^(t, z) - e iztJl Vo) - el e lzJl{t ~ 2s) (y^~ izsJl + ie izsJl J 1 V 2 ^(s, z)ds. 
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Thus we obtain 



a{t, z)iP{t, z) = e iztJ a{0, z) + e / e izJ{t - s) W s ip(s, z)d 



which yields We will show (E3D-(E6l). Using \l/a(t,z)\ ^ 2 for sup tm \V t \ < |z| and 

|V>o(MK el Im2 l', we have 

|# n (t,z)| <2\e\ f e^^-^Wt^n^z^dU 
Jo 



t r t\ 



Jo 



tn—1 



i ' 



which gives ([33]) . Estimates ([S3]) and J™ |W(t)|di < m(||V|| + ||V|| 2 ) imply Q33J and (1331) . 

We get 

#i(l,z) =£ J fs:a-Vo = e^o + £ 2 ^^o + 0(e 3 e |Imz| ). (3.9) 
Recall e t = e^* Jl . (ETIDjl implies 



K% = - e^tii^V* + Vl)e t dt = -%e x J{f - exV'(z), 



ATt^o = - / ex-tCiJi^ + 7 t ')y tet £ft = -iexJiV 3 ^) - ei / V t V t dt 



which yields (l3~7j) . Consider the second term ^ 2 = e 2 (^a _1 )Vo = + 0(e 3 e |Imz| ), 

where ^ ^ 

/ dt / ei_ t (iJxy t 2 + V/)e t _ s (iJxV s 2 + V> s cte 
'o 



20 



1 rt 

dt / e X -t 



o jo 



i 

-ex / dt 
'o Jo 



if K 2 - ie_ 2a Jx^K' - ie-v^JtV? - e_ 2t+2s KK 



' 7 T^2 



f/.S' 



which gives (13.81) . ■ 

In order to determine the asymptotics of the Lyapunov function we need the following 
modification. Substituting = —Ji/; T J into = aoif)(l, z)~ x clq we get 



^(1, Z )- L = -a Ja ^(l, z) T a 1 Ja 1 = -cJtf (1, z) 1 Jc~\ c = 7 - e*V£, e = ^~ 



T T -1 



2z' 



which yields 



L= l - (tf (1, z) + z)" 1 ) = i (*(!, z) - cJ* T (l, z) Jc- 1 



(3.10) 
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We determine the asymptotics of ty(l,z). Using (I3.6p . we have 

^ = + eKa^ipo + e 2 Ka~ 1 K a - 1 ^ + e 2 Ka~ 1 K a - 1 KoT 1 ^ + 0(e*e l lmzl ) 

= ipo + eK(I + eV + e 2 V 2 )ip + e 2 K(I + eV)K(I + eV)ip + e 2 K 3 ^ + 0(e 4 e |Im21 ) 

= ^ + ex l + e 2 x 2 + e 3 x 3 + 0(e 4 e |Im21 ), 

where 

X! = Ki[> , x 2 = (KV + K 2 )ij , x 3 = (KV 2 + K 2 V + KVK + K 3 )^. (3.11) 

We get 

2L = ^ + ^i+£ 2 ^2 + £ 3 x 3 -(/-£V 2 )J(^o + ^i + £ 2 ><2+^3) T ^(/ + £V 2 )+0(£ 4 el Im2 l) 

= (V» - + e(*i - J>c[j) + £ 2 (x 2 - Jx 2 T J + J^ T J - J^ T JF 2 ) 

+e 3 (x 3 - JxJ J + VgJxlJ - JxJjV 2 ) + 0(£ 4 e |Im21 ). 

Thus we have 

L = cos z + eL x + e 2 L 2 + e 3 L 3 + 0(e 4 e llmzl ), (3.12) 

where 

L l = Kl - J 2 KlJ , L 2 = X2 - J f 2 \ L, = \{h>- JxlJ+V 2 JxJj-JxJJV 2 ). (3.13) 

Below we need the identities 

J*V T J = -V, J*W T J = -W, J*J X J = -J x . (3.14) 
We determine the asymptotics of L. Let u t = J * V 2 ds, f t = V!,e 2s ds. 

Lemma 3.2. IfV, V G Jrff, then asymptotics (11.101) . the following identities and asymptotics 
are fulfilled: 

L 1 = Ys\nz, zeC, where V = / V 2 dt, (3.15) 

Jo 

L 2 =i^(rV> + V>y)-ismz I J^Vtdt + L 21 + L 22 , (3.16) 
^ Jo 

L 2 i = ~ (eiu'u + e_iW)dt, L 22 = ^ ^ (ei/7 + e-iff^dt, 

L 22 = o(el Imz l), (3.17) 
^(™) = H^f-^ 2 + ^i(^ + K^) + (^) 2N ), K = n™), (3-18) 



A j (^) = cos(^-^) + Q( ^" l t |n| as* = 7m + o(-Y j = l,..,N. (3.19) 
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Proof. Recall e t = e iztJl . lETIIjl . rtOjl give 

• i 



r H [ J = J*(K^ ) 1 J = J* / (ei_ t W t et) 1 Jdt = / e_ t J*WV Je^id* 
Then 



e_ t W t e t _ x dt. 



which yields (ETI511 . We determine L 2 = §(x 2 - JkJJ). Using fl3TT|) . fl4T|) we get 
r^J=rJ { e tVt T +J o {e 8 Wje a -tds)W?e X - t dtJ=j(e- t Vt+j (e. s W s e s . t ds^W t e t ^dt. 



This yields 



L 2 — F + S, S = ]^j^dtj^ (eLtWtetsWsCs + e_ s W s e s _ t W^_i) rfs 
and using VU = — i.J\V 2 — V we obtain 

F =J^e 1 . t W t V t e t + e. t V t W t e t . 1 ^ = -J^tiiJxV 2 + V^)V t e t + e_ t V t (tJ 1 V 2 + V;)e t ^ 



\ dt 
2 



exV'Vt + e-xVtV'jdt 



e_i - ei 



VjVtdt = -iJi sin ^ / VfVtdt 



since J" VtV^dt = — L V(V t dt. Consider the second term, 



dt j [ei- t 
o Jo 



J x V 2 e t . s J x V 2 + i.hV 2 e t . s V' s + %V' t e t . s J^ + Vje^X 



e s + 



+e_ 



J x V?e a . t J x V? + iJ x V 3 2 e s - t Vl + %V' s e s . t hV 2 + V' s e s - t V' t 



e t -i 



1 rt 
dt 

^0 



- e x V 2 V 2 + z ei _ 2s Ji^K + iei-2tV;j{V? + e^ 2t+2s VX 



rl\r> 



' t ir2 



+ 



e-iVX 2 + ie^nhVfVL + ie^VlJiV? + e„ 1+2i _ 2s K'K 



I J T/2 



jf dt(e 1 u'u + e^ 1 uu^dt+ l -^(yV' + V , r)+ J dtfej'f + e^ff^jdt 
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where V'{z) = V'(t)e2tdt and u = J* V 2 ds, f = f* V^e 2s ds, and here we used 
J dt J (eiV?V? + e-iV?vf)dtds = J dt^uu + e^w/jdt, 

[ (u'f + f'u + uf + fu')dt = Ul h + fmi = rv> + V"T, 

Jo 

J dt jT (ei-^tsVX + e n „ 2s _ x VX)ds = J &(e x f f + e. x ff^dt. 



We have e±\ = (—l) n l2N at z = 7m and then 



L 21 (im) = I (u'u + uu'Jdt = -(-l) n — , (3.20) 



o 



M™) = i-^- / (/'/ + //')dt = { -^f(irn) = { -^L(V') 2 (r m ), (3.21) 
which yields (13181) . Using fl332l . fl3~T5|)(l3TT8l) we obtain 

L(z) = cos(z-er)+e 2 { < L 2 (nn) + ( y -l) n —)+0( y E 3 ) as z = im + 0(l/n). (3.22) 

Recall the simple fact: Let A, 5 be matrices and and be spectra of B. If A be normal, 
then dist{a(A), a(A + B)} ^ |5| (see [Ka,p.291]). 

The normal operator cos(,z — e'f) has the eigenvalues cos(z — euj),j = 1, .., N with the 
multiplicity 2. Using the result from [Ka] and asymptotics f)3.12p and identity (13.181) we 
deduce that the eigenvalues Aj(z) of matrix L(z) satisfy the asymptotics (I3.19p . The proof 
of (11.101) is similar. ■ 



4 Proof of the main theorems 

We need the following results from [K4]. 

Lemma 4.1. Let V, V G J4f. Then the following asymptotics hold true: 

<$>(z,u) = (cosz-u) 2N + o(e Nllmzl ) as \z\ -»• oo, (4.1) 

where \u\ ^ A for some constant A > 0. Moreover, there exists an integer such that: 

i) the function 1) has exactly N(2uq + 1) roots, counted with multiplicity, in the disc 
{\z\ < 7r(2n + 1)} and for each \n\ > n , exactly 2N roots, counted with multiplicity, in the 
domain {\z — 2ixn\ < ~ }. There are no other roots. 

ii) the function —1) has exactly 2Nuq roots, counted with multiplicity, in the disc {\z\ < 
27rno} and for each |n| > no, exactly 2N roots, counted with multiplicity, in the domain 
{\z — ir(2n + 1)| < f }. There are no other roots. 
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Hi) Let in addition z/j ^ Vj for all i ^ j G uj s for some s = 1,..,N . Then the function 
p s has exactly 2N S (N S — l)n roots, counted with multiplicity, in the disc {\z\ < ir(n + |)} 
and for each \n\ > n , exactly N S (N S — 1) roots, counted with multiplicity, in the domain 
{ I z — irn I < |}. There are no other roots. Moreover, 

/sin2 + o(el Im2 l)\^(^-i) „ 
p s (2;) = cJ J , c s = [[{Vj-Vk), \z\ ^ 00. (4.2) 



All zeros of p s are given by z% ,a = {j,k),j < k,j,k G uo s andn G Z\{0}. Furthermore, 
they satisfy 

Proof of Theorem 1.2. i) We determine asymptotics (11.151) for z™ as n — > ±00, j = 
1,..,N. Lemma 13.21 yields l-z^ — 7m | < ^, as n — ► 00, j = l,2,..,2iV. Lemma 13.11 gives 
A m (z) = cos(,2 - ^) + 0(l/z 2 ),m = 1, .., iV as z = im + 0(1). For each m = 1, ..,iV there 
exists j such that Aj(z^) = (—1)™. Thus we have z^ = nn + 0(l/n). Define the lokal 
parameter \i by z = im + e/i, e = In order to improve these asymptotics of we need 
asymptotics of the ^(1, z) as z = nn + 0(l/n),n — > ±00 given by (13.71) 

z) = e izJl (l - ieJ 1 T n + 0(e 2 )) , T n = Y- iJ^, = J V'(s)e i27Tn ' hs ds 

where ^ = V 2 (t)dt = %Q)%, % = diag{^, v N }, < v x < i/ 2 < ••• < fjv- Thus we 
get 

A = [ L) ^ {l : Z} = ^ : i = J x (v n - p + O(s)) . (4.4) 

— IE —IE V / 

Hence we study the zeros of the equation 

det(^-iJ 1 V( n) + O(e)-^=0, K=(~* ^ j , < = j\\t)e-^ nt dt, (4.5) 

where /i G C. We will use the standard arguments from the perturbation theory (see 
[Ka,p.291]). Let A, B be bounded operators, A be a normal operator and a(A),a(B) be 
spectra of A, B. Then dist{a(A), a(B)} < 

Let ( m ) n ) {1) 2, •-, iV} x Z be the eigenvalues of the self-adjoint operator Y— i JiV^. 
Using the arguments from the perturbation theory (see [Ka,p.291]), we obtain that Eq. H4.5[) 
has zeros cj^, (to, n) G {1, 2, .., iV} x Z such that = + 0(n _1 ) as n — ► ±00, which 
yields f TTTT5|) . 

Consider the case ^1 < ... < i/jy We shall determine asymptotics (ll.ISp for the case 
a = (to, to), rn = 1,..,N. Let = 7m + e/x and /i — i/ m = £ — > 0. Using the simple 
transformation (unitary), i.e., changing the lines and columns, we obtain 

det (V - i Ji + O(e) - = det \^ ^ j = det A 4 det 
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At = ^ ^ _^J+0(e), b = -iv n , ta ,a = (m,m), A 2 ,A 3 = 0(6 n ), 5 n = \V^\ + \e\, 

A A = % m © r 0m + 0((J n ), % m = diag{^- - C, J ^ m}, 
X = A x - A 2 A 4 1 A 3 + 0(e) = A l + 0(0), <p=\e\ + \K\ 2 , 

which yields 

= det K = £ 2 - \b\ 2 - e&i + a& 2 + && 3 + O(0 2 ), & 1; 6 2 , , b 3 = 0(0), 

where 6i, b 2 , b^ are analytic functions of £. Rewriting the last equation in the form (£ + a) 2 = 
(c + /5) 2 + O(0 2 ), a, /3 = 0(0) and using the estimate ^ x 2 + y 2 — x ^ y for x, y ^ we get 
£ = ±c + 0(0), which yields fOHjl for i= j. 

Consider the resonances. We shall determine asymptotics fll.ISp for the case Uj ^ Uj> 
for all j 7^ j' G u s for some s = 1,..,Nq. By Lemma 13. 2^ the zeros of p s have the form 



a 



(j,f),j,f euj s ,j < f,neZ and satisfy \z^ - nn\ < tt/2. 



Asymptotics (OD]l yields A j (z)-A jl (z) = (vj -u j/ )^ + 0(z- 2 e^ m ^), \z\ -> oo. Then 

— 7m| < 7r/2 yields = irn + 0(l/n) asn-> oo. 
We have the identity Aj(z) — Aj/(z) = at z = z^. Then using (13.191) we have 

cosU ^_^_ cos ^_^U2(-irsin^^sinr^-vrn V * + ^ - 



Q 27m/ V a 27rn/ 47m V a Aim J n 2 

which yields (11.161) . i.e., 

z? = *n + e{a + + Z?), a ± = ^^ } e = JL z? = 0(8 n ). (4.6) 

We shall show that for large n in the neighborhood of each im + ea + the function (Aj(z) — 
Aj/(z)) 2 has two real zeros resonances (counted with multiplicity). Introduce the functions 

f m (ji) = 2(2nn) 2 (l - (-l) n A m (7m + efx)) = (ji- v m ) 2 + 0(5 n ). (4.7) 

For the case \i — > a + we get 

/ m (/i) = (a+ - z/ m ) 2 + o(l), m = l,..,iV, and f m (/i) = a 2 + o(l), m=j,f. (4.8) 

Hence the function /j — /j/ (maybe) has the zeros, but the functions fj — f m ,m ^ have 
not zeros in the neighborhood of the point a+. 

Note that these functions are real outside the small neighborhood of a + , otherwise for 
any complex branches there exists a complex conjugate branch, but the asymptotics (14 .7p 
show that such branches are absent. 

We have two cases: (1) let / s (aO, s = j,f be real in some small neighborhood of a + . Then 
the function fj — jy has at least one real zero since by Theorem II. II the functions fj, jy are 
strongly monotone. Thus (fj — fj') 2 has at least 2 real zeros. 

(2) Let f s (y), s = j,j' be complex in some small neighborhood of a + . Then they have at 
least two real branch points. Thus (fj — fj/) 2 has at least 2 real zeros. 
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Hence (fj — fj/) 2 has exactly two real zeros, since the number of resonances (in the 
neighborhood of the point Tin) is equal to N S (N S — 1). 

We determine the sharp asymptotics of resonances. Recall that 

(— l) n fy(l z) — I 

a = — V-^ — = -fi)- iv> + o(e), r = %@ %. 

—IE 

The operator A — a_ has the eigenvalue £o = T "."~ 1 _ a _ Q f multiplicity two, since 
r ntS = (_i)V £ ( a -+°( 1 )). The operator J^Y - a+) -a- = (% - vj) © (u f - %) has two 
eigenvalue (= 0) and other eigenvalues are not zeros. Using the simple transformation 
(unitary), i.e., changing the lines and columns, /i = a + + rGlR, r = r n>s — > 0, we obtain 

det{A - a_ - = det ( ^ ^ 2 J = det A A det K, 



~Y * l \) + 0(e), K — Ax — A 2 A^A 3 = ( ~ T . ~ * + ai ib ~ * + a * 
ib r-£ J w 4 V tb + a 3 r + a 2 



A 4 = f 0j © + 0(5 n ), b = -< >a , a = (j, /), A 2 , A 3 = 0(5 n ), 

the function ai, a 2 , 03, a 4 = 0(0), = |e| + |V^| 2 and they analytic with respect to £ in some 
small disk. The function det K has the form 

det K = £ 2 — r 2 + \b\ 2 + a x (r — £) + a 2 (— r — £) — i> va 3 ~ * a 4^ — 0,403 (4.9) 

Then = det K = (f - £ ) 2 (1 + 0(£ - &)) for £ ^ where £0 = (( ~ 1) "^ s ~ 1) - a_ is the zero 
of F of multiplicity two. Then £0 — 0(0) and we have (r — / y) 2 = (\b\ — f3) 2 + O(0 2 ) where 
7, (3 — 0(0). Then using the estimate \J x 2 + y 2 — x ^ y for x, y ^ we get r = ±|6| +0(0). 
■ 

Proof of Theorem 11.31 We consider iV + , the proof for N~ is similar. 

(i) Assume that A" + = 00. Then, due to the Lyapunov-Poincare Theorem and Theorem 
11.21 there exists a real sequence Zk — > 00 as — > 00, such that 2^ G #™j( m ) m } for each 
m = 1,..,N. Hence, ^m=i9{j( m ) m} ^ ^' Using asymptotics (11.181) and — > 00, we obtain 
v\ + = ... = un + z/j(7v). Moreover, the estimates z/i < ... < un yield z/^ > ... > z^jv), 
i.e. j(l) = N, j(2) — N — 1, ... Then, i/x + z/^v = + = which gives a contradiction, 
(ii) Let 2a = v x + v N = v 2 + v N _ x = ... Due to f|1.18|). 7rn fc + a G n^ =1 ^ +1 _ mm} as fc -> 00. 
Then the Lyapunov-Poincare Theorem yields im^ + a ^ cr(/C), fc — > 00, i.e. A^ 9 = 00. ■ 
Proof of Theorem 11.51 i) We need the following result from [K4]: Let V G ffl '. Define 
the quasimomentum kj +N = kj = arccos Aj(z) = i\ogrj(Aj(z)), j = 1,..,N, see fll.23p . 
(11.22p . Then the averaged quasimomentum k = 517 X)i % = "h Si ^ s analytic in C + 
and k : C + —* k(C + ) = WL(h) is a conformal mapping for some h G C us . Furthermore, and 
there exist branches kj,j 6 1,JV such that fll.26l) - fll.28l) hold true. 

ii) Let V,V G J$f. We need the following results from [K4]: let for some constants 
Cq, C\, 2 the following asymptotics hold 



exp —%2N ( z -9l-9l- C2 + ) ; as« = <y, y - 00. (4.10) 



det(M(^) + M" i (z)) ,. M . 

V 2 z z ^ 
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Then 

7 f \ QO Ql Q2 +0(1) ^|| r ^ n /„ „X 

= z 5 , as y > ro \x , y — ► 00, tor any r > U, (4.11) 

where Cj = Qj, j = 0, 1, 2, Q 2 = I® -\- 1% — %■ Using these results and asymptotics from 
Lemma Owe obtain ([L29])- ([1130]). ■ 



5 Example of complex resonances 

Let below N = 2. Consider the operator /C„ T = —iJ\^ + K, jT , i/ = 1, |, |, .., r G R acting in 
L 2 (R) 4 , where the real periodic potential V UtT is given by 

v " = U Y) • = - ( ri,w T6 o (t) ) • s e R + ^ N ' 6 " e C(T) ' (5 ' 1} 

We need another representation of K V)T . Recall that % = -j^(Ji + iJ) = ^* and iden- 
tities (12.221) give JC~ T = ^JC VtT & = jA — V VjT . Using the unitary transformation y = 

(j/i, 2/2, 2/3, 2/4 ) T ^ Uy = (2/1, 2/3, 2/2, 2/4 ) T in L 2 (R) 4 we define the new operator P U>T = 
A 



UlC-^r = U+W u , T , where 



w„ = - ^ = ( * 2 rt - ) ■ t = ( ! J )■ J=(_i J ) • <"> 

We rewrite W„ )T in the form W T)V = W° + rb u J 2 , where W° = a jj ^ . If r = 0, then 

we have the unperturbed operator P° = j'4 + with a constant potential W°. 

1. The 2x2 Dirac operator. We consider the simple example of the 2x2 Dirac 
operator P® = jj| + aj 2 , a > acting in L 2 (R) © L 2 (R). The spectrum of P° is purely 
absolutely continuous and consists of two intervals (— 00, — a), (a, 00) separated by the gap 
(—a, a). In this case only one gap is open and other gaps are closed. The solution of the 
system jy' + aj 2 j/ = zy has the form y = e ±tkt yo, for some constant vector y G C 2 and k 
satisfies 

det(i/cj + aj 2 — zl 2 ) = —k 2 — a 2 + z 2 , fc = = \J z 2 — a 2 , 

where the quasimomentum fc : C \ [—a, a] — ► C \ [ia, —ia] is the conformal mapping with 
asymptotics k(z) = z + o(l) as |z| — > 00. Note that 

A = j(z - aj 2 ), A 2 = - aj 2 )) 2 = -z 2 + a 2 - za}j 2 ~ zaj} 2 = a 2 - z 2 = -k 2 . (5.3) 

Then the fundamental solution of the equation + aj 2 vp® = zip®, ipx(0, z) = h is given by 

z) e~* A T { - tAr T( itA)2n {tA)2n+1 ) 
W,z)-e -2^ , ZA (2n) , (2n+ iW 
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^/(-l) n (tk) 2n A(-l) n (tk) 2n+1 \ , A. , 

> - — - v / v '— — = costk- -smth. 5.4 

^\ (2n)\ k 2n + l / k 



n>0 



Thus the Lyapunov function has the form 



Tie 



-A 



A°(z) = = cosfc(z), fc = fc(z) = V* 2 - a 2 , [-a, a]. (5.5) 

2. The 4x4 unperturbed operator. Consider the 4x4 operator P° = j-| + PU in 

L 2 (M) 2 © L 2 (IR) 2 , where P^ — f °q 2 g ^ is tri e 4x4 matrix. We rewrite this operator 

in the form P° = P® © P®, where the Dirac operators P® = j-i + a] 2 and P 2 ° = j4 act 
in L 2 (R) © L 2 {R). The corresponding equation j-^ip + W°ifj° = zip has the fundamental 
solution il}°(t) given by 

^(t)=^(t)(B^ 2 (t), ^ = e~ A \ iP°2 = e-i zt , A = j(z-aj 2 ), t>0. (5.6) 

Thus using (Q ]) - (j5TBjh ffTTTTD we obtain that D°{t, ■) = det(^°(l) - rl 4 ) satisfies 

D°(r, ■) = (r 2 - 2A?r + l)(r 2 - 2A°r + 1), A^ = (T° - (-l) m y^)/2, 

T° (z) = Tr(e~ mA + e- j " m ) = 2(cosmA;(z) + cosmz), p°(z) = (cos k{z) - cosz) 2 , (5.7) 

for any z G C, to = 1, 2. We will determine the zeros of p°(z). We have = cos A; — cos 2 = 
2 sin sin k = k(z). Then we obtain k± z = 27m, n6Z \ {0}, which gives zeros r° 

(each zero has multiplicity 2) of p° by 



r£ = 7m+A fc(r°) = | 7rn * |n | > 2 } , nGZ\{0}. (5i 

4 ™ -™ + £- if n < £ 



We determine the periodic spectrum for the equation jy' + W°y = zy. Using (I5.6p . (l5.7p we 
have det( , o (l, z)^I^) = 4(cosfc(z)=Fl)(cos2:=Fl), which yields the periodic and anti-periodic 
spectrum multiplicity 2 



/ a 2 

z 3 ' 1 = to, n 6 Z, and z\ = irn\ 1 H — -, nGZ \ {0}, (5.9) 

V 7i n 

and Zo,2 = ^ a nas multiplicity one. Note that the zeros z° 7^ for all n, m G Z, p = 1, 2. 
In this case there are no gaps in the spectrum. 

3. The perturbed case. Consider the 4 x 4 operator P TiU =jf t + W T ,v in L 2 (1R) 2 ©L 2 (R) 2 , 
where the 4x4 potential W TjU = W° + rb u j 2 J 2 satisfies (15 .11) . (jl.20p . We show that there exist 
the non-degenerated resonance gaps for some W T}U . In this example some resonances are real 
and some are complex. The fundamental solution ip T,u (t,z) of the Eq. jip' + W TiU ip = zip 
satisfies the integral equation 

'.0/+ „\ 1 „ / „/,0/ 



^"(t, z) = ^(t, z)+r I ^ U (t - s, z)iK(s)j 2 J 2 ^(s, z)ds. (5.10) 
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Then ip T,,/ (t,z) has asymptotics 

tp T ' u {m,z)=^ {m,z)+Ttfj 1 {m,z,u) + T 2 ^ 2 {m,z,u)+O{T 3 e mllrazl ) as r -> 0, (5.11) 
uniformly in v = 1 A .., m = 1,2, z G C, where 

rl)\m,z,u)= ^{m-t,z)b v {t)^J 2 ^{t,z)dt, (5.12) 



ip(m,z,v)= / i/j°(m -t,z)b u (t)} 1 J 2 dt ip°(t - s, z)b iy (s)i 1 J 2 ip (s } z)ds. (5.13) 
Jo Jo 

The identity Tr?/> (t, z)\\J 2 = 0, (i, z) 6 E X C yields Tr ^(t, z, z/) = 0. Thus we obtain 

T£"(z) = Tr ^(m, z) = T° m (z) + r 2 T m2 (z, z/) + 0(r 3 e m l Im ^), m = 1, 2, (5.14) 



t 

0/ 



T m2 (z,v) = b u {t)dt I b v (s)F m (t,s,z)ds, F m (t, s, z) = Tr ji J 2 vp (y, ^)ji J 2 ^o(C, ^) ? 
Jo Jo 

(5.15) 

uniformly ini/ = l,|,..,|r| < 1, z G C, where y = t — s,C = m — y. 

Lemma 5.1. Lei V^. v satisfy (11.191) . (11.201) . Then the following asymptotics hold true 

T 12 (z,u) = ?M + ( c l*M )j T 22 (^,z/) =T 2 °(z)+40(z) + o(e 2 l Imz l), (5.16) 

T^) = + + o(r 2 el Im ^l), (5.17) 

T^z) = (1 + r 2 )T 2 °(z) + r 2 40(z) + o(r 2 e 2|Im ^), (5.18) 
p^( 2 ) = - Tll^L = (1 + r 2 )p°(z) + r 2 2(^) - 1) + o(r 2 e 2 l Im2 l), (5.19) 



det(^ T '"(l, =F h) = (l + r 2 )D°(±M)+T 2 (r°-l- ^ 1 + o(e 2|Im ^)) , (5.20) 
as i/ — > ; uniformly on |r| <1,zgC, and where <fi(z) = cos /c(z) cos z + sin 2; sin k(z). 
Proof. Using b v — > 5 per = J25(t — ^ — n) and •) G C(M 2 ), j = 1, 2, we obtain 



Jo Jo 



K(t)b u (s)F 1 (t, s)ds = iFx (± i) + o(el Im ^l), (5.21) 

J'dtj\ v (t)b v (a)F 2 (t,s)ds = ~^(^~) +iF 2 (|,2) +F a (2,l) + (e»l^l) (5.22) 
as z/ — > 0. Thus, if m = 1, then (15.211) gives 



2i 2 (z,^) 



i Tr j x J 2 ^°(0, z)j! J 2 ^o(l, *) + o(e |Imz| ) = ^ + ( e l i-l), 
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If m = 2, then F 2 (§, |) = F 2 (§, |) = Tr^o(2, *) and the identity J 2 ^° = ( ) y ields 



V>" o 

'4 s) = ^ ( tftofS > ( «<0 1° ) ^ 

A A 
= 2 Tr ji(cos z — j sin 2)^ (cos — — sin k) = 2 Tr(cos;z + j sin 2) (cos — — sin fc) 

A 2! 

= 4 cos z cos k — 2 Tr j — sin z sin = 4(cos z cos + — sin z sin fc) = 40. 

2fc fc 

Then (ET22p gives 

T 22 (^,z/) =Tr^ o (2,^)+40(z) =T 2 o (z)+40(^) + o(e 2|Imz| ) as 1/ -»■ 0. 
Substituting fl5TT6|) into (15TT4D we obtain (OTjl . fl5TT8|) . The asymptotics floTTTI) . (I5TT8D imply 

p ^( Z ) = ^Og) + 4 - ^1^)! = (1 + r 2 )p (z) + r 2 2(0(z) _ 1} + o(r 2 e 2|Im,| ); 

and D T ' U (±1, •) = det(^ T '"(l, •) =F h) =(T 1 T ' I/ - l) 2 - p r > y satisfies 

D^(l,z) =((1 + - I) - (1 + r 2 )p°(^) - ^(0(z) - 1) + o(r 2 e^) 

= z) + r 2 - 1) - p°(*) - {<P{Z) ~ 1] + o(l)) 

as v — > 0, which gives (I5.20p . The proof for D T,U (— 1, •) is similar. ■ 
Proof of Proposition 11.41 We have the simple asymptotics 

2 2 
k(z) = z- — + 0(z- 3 ), cosk(z) = cosz + — smz + 0(z~ 2 e llmz[ ), (5.23) 

p°(z) = (cos£;(z) - cosz) 2 =(^^-^ + 0(z~ 3 e 2 l Im ^) (5.24) 

as \z\ — > 00. We take x > such that the disks B n = {\z — r°| < x},n G Z are not 
overlapping. Define a constants rj n = min| z _ r o| =>e . |p°(z)| > 0. Thus using (I5.20p we get 

\p T ^(z) -p°(z)\ = r 2 0(l) = — p (z)O(l), \z-r° n \ = xas r -> 0, (5.25) 

for each n. We also obtain for |z| = 7rn + 1 the following estimates 

\p T ' u (z) - p°(z)\ = r 2 e 2|Im ^O(l) = r 2 |4sin^| 2 0(l) = T 2 \z\ 2 \p°(z)\0(l), n -> 00. (5.26) 

Thus we take large n and sufficiently small r such that |p(z) — p°(z)\ ^ ||p°(z)| on all 
contours \z\ — it (no + 1) and \z — r°| = x, |n| ^ no- Then by the Rouche theorem, p has as 
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many roots, counted with multiplicity, as p° in the disks {\z\ < min + 1} , {\z — r®\ < x}. 
Since p° has exactly one double root at r°,n 7^ 0, and since n can be chosen arbitrarily 
large, we deduce that in each disk {\z — r°| < x}, 1 ^ \n\ ^ no there exist two zeros r^ u (r) 
of p T ' u for sufficiently small r, v. 

Consider the zeros v^ v {r) of p T,u in the disk {\z — r°| < x} for fixed n, 1 ^ n ^ no- The 
proof for the case n < is similar. Recall <p(z) = cos k(z) cos z + sin z sin k(z) . Consider 

<p(z) — 1 at the point r° n = nn + x n ,n ^ 1, where x n = Using (15. 8p we obtain 
0(r°)-l = cos 2 x n -l + -^-sin 2 x n = f-^r--l) sin 2 x n > 0, if 1 ^ n < (5.27) 

0(r°)-l = cos 2 x n -l--^sin 2 x n = -(l + -^)sin 2 x n <O, i/ n > (5.28) 
We rewrite the function p T)V in the disk {|z — r°| < x} in the form 

^r) = ^^ = (,-r°) 2 /^)+r 2 1 ( 2 ,r), /(*) = z 6 S„, (5.29) 

for sufficiently small fixed u, r. The functions R(z,t), f(z),(f>i(z,T) are analytic in (z,t) G 
-Bn x {|r| < e} for some small x, e > and satisfy 

/(r°)>0, 1 (r°,O) = 2(0(r°)-l)+o(l) as !/ - 0. (5.30) 

Applying the Implicit Function Theorem to R(z,t) = and using fl5.27p - fl5.30p we obtain 
a unique solution r^ u (r) of the equation $(r^(r),r) = 0, r G (— r , r ), v (0) = r°, for 
some r > and here r^^r) is an analytic function in {|£| < r } and satisfies fll.2ip with 
R n = — 2(0(r°) — l)//(r°). The proof of the statement ii) is similar. ■ 

6 Appendix 

Lemma 6.1. Let V, V G J$? and let r > 0. Then fory ^ r|x|, y — > 00 following asymptotics 
hold: 

TrLi(z) = ||\/|| 2 sinz, (6.1) 

/ Tr^ 2 \\V'\\ 2 + o(l)\ 

TtL 2 (z) ={iH, - + "J 1 J ) cosz, (6.2) 

TryL 2 (z) = cos^Tr^d + G - —J + o(e |Im2: |), (6.3) 



TrL s (z) = icosz 



Ti:(jf v^dt + Gi + G-^j + o(e |Imz |), 

G 1 = 1 / [ J^Vtdt, G = V I J lU ' t u t dt, u t = [ V s 
Jo Jo Jo 



(6.4) 
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detiW = 2 -e X p-^-£-^-^), (6,) 

where H = Tr £ V t 2 dt, Hi = Tr -U^Vtdt, H 2 = Tr (V/ 2 + U t 4 ) dt. 
Proof. Identity (13.151) implies (16.11) . We will show (16. 2p . Asymptotics ( 13. 12[) yields 

L 2 = -isinz [ J^Vtdt + i^ifV' + uV)+L 2 i + L 22 , (6.6) 
Jo * 



LlOA — 



as \z\ — > oo,y ^ r|x|. We determine Tr L 2 . Using (I2.12j) we have TreiJi^V + V 7 "^) = 
and then 

Tr L 2 = -i sin zHi + Tr(L 21 + L 22 ). (6.7) 

Using (I2.15P we get 

TrL 21 = -y^ 1 (e lM 'u + e_ 1 w)rft = Tr £ (u'u + nu^di = - C ^Trf 2 . (6.8) 

Due to (I2~T5|) . (12T2|) we obtain 

TrL 22 = — [ (e l ff + e„ l ff)dt = Tr [ [ e izJl{1 - 2t+2s) V^dtds (6.9) 
2 Jo ^ ' Jo Jo 

and using 

* 1 f jMt-s) hrsdtds = f 1 if&M + £(i) ; Z" 1 f t e -^ t -s) ftredtda = 2^!^l f (6 .io) 

) Js Jo 2z z Jq J s z 




for f,h& L 2 (0,1) and Tr Ji(U') 2 = 0, we have TrL 22 = ^(jl^T + o(l)Y which give 
(16.21) . We will determine Tr^L 2 . Using (I6.6P we have 

Tr fL 2 = Try(^ — isinzJ^ J x V'Vdt + L 21 ) + o(e |Im " 1 ), (6.11) 

since Tr fei Ji^U 7 + vV) = and L 22 (z) = o(e\ lmz \). Due to e izJl = cosz + i^sinz and 
$l{u'u + W)cft = r 2 , Tr r 3 = we get 

m j ? t cos 2; - sin 2, T Z" 1 . , ,, , cos,?: ., q _ 
Tr rL 21 = Tr r 3 - i Tr r Jj / (it'u - W)dt = Tr r 3 - z sin zG 

2 2 y v ; 2 

and together with (16.111) we obtain (16.31) . since sin z = zcos+0(e~' Im2 '). 
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We will determine TrL 3 . Using (13.131) we rewrite TrL 3 in the form 

4 

Tr L 3 = Tr x 3 = Tr(KV 2 + K 2 V + KVK + K 3 )ip = ^ A k , 

1 

A x = Tr KV 2 ifj , A 2 = Tr K 2 Vip , A 3 = Tr KVKtfj , A A = Tr K 3 ^ 0} (6.12) 

and recall (Kf)(t) = f* e t _ s W s f(s)ds } W = -iJ x V 2 - V . Using fF2~T2l . fT2~TD . (j2TT5|) and 

e x = cos z + iJi sin z we have 

A x = Tr I e x W t V 2 dt = Tr [ e x (-iJ x V 2 - Vl)V 2 dt = Tr [ -U x e x V 4 dt = 



= Tr — iJ x (cosz + iJ x sin z)V t 4 dt = sinz Tr / V t 4 dt. 
Jo Jo 

The similar arguments give 

A 2 = Tr [ [ a-t+sWtetsWsVsdtds = Tr [ [ e^ t+s (iJ x V 2 + Vl)e^ a {iJ x V? + V^)V s dtds 
JoJo JoJo 

= Tr jj U^VfVX - e^ 2t+2s VX 3 )dtds (6.13) 

and 



A 3 = Tr/ / e 1 -. t+s W t V i e t - a W a dtds = Tr [ [ e X - t + s (iJ x V? + V^V t e t ^(iJ x V 2 + V' s )dtds 



'OJO JOJO 
-I nt 



= Tr J] iJ x (e x V;V t V 2 + e x „ 2t+2s V 3 V^dtds. (6.14) 
Summing (I6T31) . (101)1 we get A 2 + A 3 = F + F u where 

F Q = Trf [ U x e x (v 2 V' s V s + VlV t vf) dtds = Tr U x e x rf V^Vtdt = icosz Tr J x y[ VjV t dt, 
JoJo^ ' Jo Jo 

(6.15) 

since iJ x e x = iJ x cos z — sin z and Tr V J VjV t dt = 0. We will show 



1 r l 



F x = Tr I I *J iei _ 2t+2s ( - U/K 3 + V t 3 V:)dtds = o(e |Imz| ). (6.16) 



We use the standard arguments. If V, V, V" G Jff, then integration by parts gives f 1 6 . X 6 1) . If 
V, V G Jf, then there exists P h , P' h , P£ G such that \\V - P h \\ + \\V - P' h \\ = h for some 
small h^O. Then F x = o(e |Im2| )(l + 0(h)), which yields flfU6]) . since h is arbitrary small. 
The similar arguments give 



ms 
e x ^ t+p W t e t _ s W s e s _ p W p dtdsdp = F 2 + F 3 , 
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where 

»i r t r p 



F 2 = Tr f [ f (UJeiVfVfVfdtdsdp, 
Jo Jo Jo 

F 3 = Tr J J* J\iJi) (e l+p „ s V t 2 VX + e x ^ 2t+2p VXXv + ei-^VX^) dtdsd P 
Using u s = f* V t 2 dt and V t 2 dt = Y — u s we obtain 



f f f sinz 

F 2 = Tr(i Jx e i) / / u' t u' s u s dtds = Tr(iJiei) / i^f — u s )u' s u s ds = i cos z Tr G — Tr Y 3 , 

JoJs Jo 6 

(6.17) 



since 



Tr(iJ iei ) t u s u' s u s ds = Tr + UsU ' s u s + u X)ds = Tr ^^ 3 = Tr T 3 

Jo 3 Jo 6 6 



and 



Tr(iJiei)'f / u s u s ds = i cos z Tr / u' s u s ds — sin zTr ^ / u s u s ds 
Jo Jo Jo 

sin £ /* sin £ 

= % cos 2; Tr G — Tr Y I (u' s u s + u s u s )ds = i cos z Tr G — Tr y 3 

2 j 2 



where we used: Tr J\V 3 = and Tr ABC = TtACB for real self-adjoint matrix and real 
representations ( 12.231) . (12.261) of V. Using standard arguments (see the proof of (16.161) ) we 
have F 3 = o(e' Im2 '). Summing Ax, ..,A A we have (16 .4p . 
We will determine (16.51) . Asymptotics (13.121) yields 

' I 2N + S, S = iEf + — + ^ + 0(e A ), S = 0(e) (6.18) 



cos z COS z COS z 

as Imz — > 00. In order to use the identity 

det(/ + S)=e*, $ = TrS-Tr^ + Tr^ + 0(e 4 ), \S\ = 0(e), 

we need the traces of S m ,m = 1,2,6. Due to (IBTTgj) . we get Tr ^ = -ie 3 Tr f + 0(e 4 
Using fl6TT8|) . fl6TT]) -(l6T4l) we get 



_ Tr — = Tr e 2 — - ze 3 ± + o(e 3 ) ) = Trie — ie iGi+G — + o(e 3 ) , (6.19) 

2 V 2 cos z / \ 2 V 2 / / 

/ Y 2 \ / Y 3 \ 

TrS = ieH + e 2 Tr [iHi - e 2 — J + ie 3 (H 2 + G 1 + G - — + o(l) J , (6.20) 

and summing (??)- (l6.20j) we get $ = ie7i + ie 2 ?^i + ie 3 7i 2 + o(e 3 ), which yields (16. 5p . ■ 
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